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ABSTRACT 


eT expected-cost solutions to the single-period 
inventory problem (the Newsboy problem) have been well known 
for many years, for risk cases where the distribution of 
demand is known. Also well known are minimax cost and 
Laplace solutions for the uncertainty case where only the 
range of demand is known. This study explores the use of 
order-statistic-based quantile estimators as decision pro- 
cedures while data gathers during the early decision periods. 
Simulation, using a variety of demand distributions and 
unit-cost values, provides recommendations on which period 
to leave the minimax procedure for a suggested data-based 
decision rule which will be asymtotic to the optimal risk 


procedure. 
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I. INTRODUCTION 


The classic single period inventory problem is called by 
various well-known names such as "Newsboy problem," '"News- 
paperboy problem," and so on. [6]. In this thesis it will 
be called the Newsboy problem. 

Solutions to the Newsboy problem are well known under 
both "risk" and "uncertainty" conditions. [7]. The purpose 
of this thesis is to study the Newsboy problem under initial 
conditions of "Uncertainty," (1i.e., when the possible future 
for demand is not yet known) and then, for the subsequent 
periods as data gathers, to suggest, structure, and evaluate 
decision procedures during this transition phase to decision 
under "Risk." 

The inventory model studied in this thesis is a time 
independent (stationary) and linear-cost type. The study 
will consider the problem of minimizing the expected cost. 
Since with the Newsboy problem both the expected profit 
maximization and expected cost minimization cases have 
identical solutions, the study should be applicable to both 
interpretations of the problem. 

Chapter II reviews the Newsboy problem for decisions 
under risk and decisions under uncertainty, and then suggests 
possible alternative procedures for "from uncertainty to 
risk" conditions. Among these, Order Statistics (a non- 


parametric approach) appear promising and are discussed in 





detail in Chapter III. Chapter IV evaluates the alternative 
nonparametric transition phase procedures of the Newsboy 
problem and as a result proposes a transition phase pro- 
cedure. In Chapter V conclusions and recommendations for 


further work are given. 





II. SOLUTION OF THE SINGLE-PERIOD INVENTORY (NEWSBOY) 
PROBLEM UNDER RISK AND UNCERTAINTY CONDITIONS 


In this section the Newsboy problem is examined under 
risk and uncertainty, and optimal decision rules for these 
conditions are given. These two cases form the bounds of 
the transition phase and permit us to characterize the 
transition phase decision problem we intend to study. 

Here, when it is said, "under risk," it is meant that 
the probability distribution of demand is known. When the 
distribution of demand is not known, and data in hand is 
not sufficient to estimate the possible futures of demand, 


then the problem is called "under uncertainty." [7]. 


A. SOLUTION OF THE NEWSBOY PROBLEM UNDER RISK 

Let Ey) be the probability density function of the 
demand for continuous populations, and Pp x) be the proba- 
bility mass function for discrete populations. If oe 1s 
the stock-out or shortage cost per unit of unsatisfied 
demand, and C. the carrying or surplus cost per unit of 
excess supply, then the cost for a particular period of 


Newsboy problem will be: 


C.(S-X), 1f 0 < X < S, and 
Cost Gls) 4c ex), if X> S, 


where X is the realization of demand for that period, and 
S is the inventory level in hand. Here S represents the 
decision variable for this problem. 


y 





The expected cost E[C(S)] which results from keeping 


an inventory level S in hand, will be: 


©O 


S 
E[C(S)] = 2 Co(S-X)pp(X) + 2. Coo(K-S) Pp) (1) 


for discrete demand distributions and, 


8 co 


EGS) 1 = Ce 1 (S-x} £(x)dx + Cae (x-S) fp(x)dx, (2) 
0 
S 


for continuous demand distributions. 

Now, the expected cost for a period is to be minimized 
and we wish to find an optimal value of S, say So» that will 
yield the minimum expected cost. In other words, (1) and 
(2) are to be minimized with respect to S. 

In the discrete demand distribution case, two necessary 


conditions for a minimum at So are: 


E[C(S,+1)] - B[C(S,}] > 0, and 


V 
em) 
® 


E[C(S,-1)]} - E[Cc(S,)] 
These conditions yield the well known result [7]: 

oes) as — : S 3 
Pp(S,-D) s ceeieos S Pp lS.) > (3) 


as conditions for a minimum expected-cost solution at So: 
Similarly, for continuous demand distributions, to 


minimize (2), the first derivative of E[C(S)] is taken with 


10 


“oe 





respect to S, and after it is set equal to zero, yields 


C 
, SO 
F,(S,) = Go @.. | (4) 
where FL (S,) is the cumulative distribution function (CDF) 

of demand at So: Thus, the 7 values which satisfy equations 
(3) and (4) are the optimum inventory levels to keep on hand 
in order to minimize the expected cost. 

This solution (under risk) is quite idealistic, since it 
assumes that the distribution of demand is known. In fact, 
in reality we are seldom certain about the distribution. In 
most cases, neither the distribution of demand, nor an ade- 
quate amount of data to provide a good basis for estimation 
of the distribution is initially available. The situation 
when: 

1. Distribution of demand is not known, and 

2. Enough data from the system to estimate the distri- 
bution of demand is not available, but 

3. We do have some idea about the range of demand D, 


say 0 < D< Dee is one we shall call "uncertainty." 


B. SOLUTION OF NEWSBOY PROBLEM UNDER UNCERTAINTY 

It has been shown that whether minimax cost, minimax 
regret, or the Laplace approach, assuming a uniform distri- 
bution of demand is chosen as a principle of choice, the 
optimal decision under uncertainty is the same [6], [7], 


namely, 


Jul 





C 50 Dax 


So = CeCe oe for continuous demand, 
Cc SO 
Cc D +1) 
~ < £s0maxt) 2. ,, 
O C. + C O : 
€ so 


for a discrete demand distribution. The results of these 
uncertainty procedures are not affected by newly generated 


demand data unless it causes one to revise the range of demand. 


C. TRANSITION PHASE OF THE NEWSBOY PROBLEM 

Now the two bounds or limit solutions (uncertainty and 
risk conditions) for the transition phase of the Newsboy 
problem are reviewed. Implicitly, to proceed in accordance 
with these current solutions, one begins under Uncertainty 
and uses the Uncertainty optimum inventory level So (based 
on the range of demand), until demand data is adequate or 
complete to estimate the distribution of demand and change 
to the optimum inventory level 5 for the Risk conditions. 
As mentioned before, the purpose of this study is to find a 
good transition phase procedure for the intermediate periods 
when some data has accumulated and is available. A transition 
phase procedure may be characterized by two decisions: 

1. When to begin the transition phase (i.e., at what 
period should we leave Uncertainty procedures to begin to 
use an estimator of So): 

2. What estimator of So to use. 

After explaining the nature of the transition phase in 


the next chapter, three candidate estimators for the optimum 


12 





inventory level So to use during transition phase will be 
proposed. 
As we remember, the optimal choice of S under Risk was 


based on a quantile. Let us call, 


_ SO 
De Cree)? 


and then from equation (4) we have, 


- lige rede 
So oa Ky - Fy (p) 5) 


which provides a minimum expected-cost solution under Risk. 
Thus, the problem is to locate the pon quantile of the 
Cumulative Distribution Function, F(X). As more and more 
data is collected, one might be able to estimate, 
g,-%,- Fle, 

as a function of the data. Some alternative approaches to 
estimate this quantile were considered to be: 

1. To consider similar systems for which we have demand 
data, and then to apply a Bayesian approach. 

2. To decide about the form of the distribution of 
demand and carry out a maximum likelihood procedure to find 
the parameters. 


3. To apply a nonparametric approach. 


For the first two cases, if one is able to decide about 


the form of distribution such as lognormal, gamma, compound- 


h 


poisson process, beta, etc., then the po quantile of demand 


13 





might be estimated quite well subject to data. It is the 


third approach, that of nonparametric statistics, which is 
the center of interest in this thesis, and forms the basis 


for the chapter that follows. 
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III. A NONPARAMETRIC APPROACH FOR THE NEWSBOY 
PROBLEM, TRANSITION PHASE 
Cimeetteass ce tOnmchc apOssl Dll I ties tor applying a non- 
parametric or distribution-free approach to the problem of 
the transition phase of the Newsboy problem will be explored. 
After explaining the reasoning for a nonparametric approach, 
three candidate estimators of the optimum inventory level 


Ee (or x will be proposed. 


A. WHY A NONPARAMETRIC APPROACH? 

As mentioned at the end of the preceding section, if the 
decision maker is able to predict the form of the demand 
distribution, then a parametric approach can be more efficient 
and have a higher convergence rate to the true value of the 
estimator than a nonparametric equivalent. Hadley and 
Whitin [5], give various considerations relating to choice 
of a demand distribution and about inventory demand predic- 
tions, in detail. For a parametric approach, such as would 
Pesult from fitting a theoretical distribution to data, they 
Suggest to use poisson distribution when the lead time is 
relatively low, or normal distribution when it is not rela- 
tively low, etc. 

Since it is quite difficult to specify the form of the 
demand distribution, a nonparametric approach was considered 
to be of particular interest for this problem. Also, Hadley 


and Whitin point out that, for single-period inventory models, 


ile) 





like the Newsboy problem, the empirical distribution function 
(nonparametric) is reasonably applicable, unless we are forced 
to operate the model at the tail of the distribution. OAs 

Performance of the model decision rules was assumed to 
| be critical during, say, the first 50 periods when data is 
inadequate. The decision maker needs to operate the system 
well during that critical time. Later on he will have a 
sufficient amount of data to treat the problem as one under 
Risk and simply use the optimal decision rules for Risk. For 
these early periods (transition phase) a nonparametric approach 
was used. 

B. THE EFFECT OF THE ESTIMATORS OF OPTIMUM 

INVENTORY LEVEL ON EXPECTED COST 

Before we suggest some estimators of optimal inventory 
level, xy? it may be useful to see how those estimators 
affect the total expected cost in a period. In the work 
that follows, it is assumed that the demand process is time 
independent (stationary), which means that seasonal effects 
and trends are removed. Also, we assume that for the previous 
period we know what the demand was. 

Before we continue on with the mathematical development, 
it is useful to expand on the conditions for an optimal 
expected-cost solution. When the demand distribution is 
assumed known, we remember from equation (2) that the expected 
cost for the case of continuous demand is: 


S © 
EL Gls) = C. f (S$ -x) fy (x)dx + Co, : (x - S) £, (x) dx. 
O 
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To find an optimum inventory level So? taking 


E*[C(S)] - SECC) iets (5) 
O 


and setting it equal to zero (necessary conditions for 


Minimization), we obtain: 


Cee Caner e es Cog | FCS = 00, > 08 
ail 
PaaS = = po 6Cand «| SS. = Far (p). 
Do Cet OS : 
Similarly, 
e~tc(s)} = ECS) = c) £4¢8,) +c. “£5 (8,) 
4 §2 SO D0 C D*"o 
S=S 
O 
(6) 
= (G2n07 Ga) £,(S,) Pac 0; 


which forms a sufficient condition for minimization. Here, 


h 


So is the value that corresponds to the pt quantile of the 


population, say Be 
We are interested in the case where the population is 

not known, but there exists a sample of size n of demand 

data, acquired during the last n periods. 

| X 


etek greeees x, be realization of demand during the 


ys 
last n periods. Although we don't know the value of x yet, 
we do have a chance to estimate this value as a function of 


available data, say: 


os _“*-] ‘nalts 
x, a G(X) .X5,----sX) @ Fo (p) . 36 (7) 
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To express the expected cost as the function of x let 


ae, 7 7, * % > 5, 


Equation (2) then becomes: 


+ 6 


A S 
E[C(S,)] = C it Cree) fy (x)dx + C2, 
0 


(x= Sy - 8) Ey (x) dx (8) 
O 


+ 
oe) 
“A-~ 8 


If equation (8) is expanded using the Taylor series around 


S it gives 


(a is 


A 2 
E(C(S,)] = E[C(S,)] + B°(C(S,)] E{é] + E~[c(s,)] ELL 


Since E“[C(S,)] = 0 by equation (5), and having E“[C(S,)] 

from equation (6), then for three terms the Taylor expansion 
oS 3 : 
om) Ea(S)) E[S"] (9) 


SO C fe 


We 


E[C(S,)] = E[c(S,)] + (C 


Equation (9) shows that: 
1. Under Uncertainty conditions the expected cost for 
a period is linear in the mean square error of the pe 


quantile estimators, 1.e€., in 


Z ae 
E[sé = B/(X - X 
[8°] = EL(X, - X,)°] 
2. If one is able to find a2 consistent estimator. 


X(m), of xy? then when n gets large the mean square error 
will converge to zero and E(C(S,)] will converge to E([C(S.)], 


providing that the moments of be exist and behave properly. 
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3. The optimal performance for an estimator here is to 
yield the expected total cost under Risk, E{C(S,)]. 

Now, we wish to find a good consistent point estimator 
x, =O Er Xn Also, it is not required, but it is Aeeaabie to 
obtain one that is an unbiased estimator of Xn if so, 
E{ (x, - x) 71 will be equal to the variance of estimator x 
Then, by equation (9), if it 1s possible to improve this 
estimator reducing its variance, it will directly reflect 
and will decrease the expected cost accrued by using that 


estimator. 


To obtain candidate nonparametric estimators of X 


? 


p 
order statistics will be used. In the following section 


we shall develop three eer extimators of Xn using 

an order statistics approach. Additional details about order 
statistics and empirical distribution functions related to 
quantile estimation are given in Appendix A. 


1. The First Candidate Estimator of 
Optimum Inventory Level: Xr) 


Let X < X ape ak be a set of n reali- 
(1) (2) (n)’ 
zations from an unknown distribution in order statistics 


th h 


form, where Xr) ES. the: +r largest value (r° order 


Statistic). Let r be defined in the following way: 
np , when np is integer 
[np +1],. otherwise, 
where p is the quantile, n is the number of periods so far, 


and [X] denotes the integer part of X. For large n, Xr) 


1S a consistant and unbiased point estimator of x (Appendix A). 
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In Figure 1, the rth 


order statistic Xr) is represented 
graphically with the distribution function F(X) and the 
empirical distribution function Gy CX). 

Know results of this estimator include the following. 


For a uniform (0,1) distribution, 


15 


E{X (5! = ELU yy] ~ hel? . . (10) 
and 
cov{[X x ] = cov[U U ] = r(n-s+1) fOr rs. 's 
(r)’" (s) Gs = Cs) Gin a 2), 
(ise) 
Similarly, for any other distribution, 
E(X-.] ¥ Fp? (5) (12) 
(Ge), Maks < 
and 
One | rinsed), a, 
(rye Gs)" (nei) 4 (42) £) (Fp (27) + fy (Fp GS) 
(15) 


provided r< os. 
Later these results will be used to compare the 


efficiency of the estimators of x, 


2. The Second Candidate HHS SCORE Ane SSS 8 SS mete 
of Optimum Inventory Level: X Optimum Inventory Level: Xr) 


As a second candidate estimator of X > Xr) is 
offered to smooth the empirical distribution function by a 
linear interpolation. It might be considered that the area 
between Gp (X) and the actual distribution function FO) 
may be a measure of deviation. So, here interpolation was 


done at the midrange, yielding a new empirical distribution 


20 





- ¢x)%9 uotzoung uotingtaistq Teotudug 


pue (xy Fy uoT}JOUNT UOTINGTAASTG sATZeETNuUND 


(T+4),, =) (T-4) 


(X* pueusq) 


X 





eect Urs )ikal 


( *eTT}UeNdH) 


SPORE! 


os) 


7a) 





function R, which would possibly have a smaller area between 

R and F(X), than the area between Gy (xX) and F,)(X). A possible 
realization for X ery is shown in Figure 2. Since the area 
between those distribution functions is subject to change 
because of randomness, the results for that case would be 
evaluated numerically rather than analytically. 


3. The Third Candidate Estimator _ 
of Optimum Inventory Level: Xr) 


As a third candidate estimator of Xo a pooled order 
statistic, combining other order statistics, 1S proposed. 
For such an approach to be successful by equation (9), it 
has to reduce the mean square error to a lower value than 


the rth 


order statistic does, if it is possible to make this 
pooled estimator an unbiased estimator of Xn Then the 
problem is one of comparing the variance of this estimator 


th order statistic. 


with the variance of the r 

To examine this subject, it is reasonable to begin 
with the uniform (0,1) distribution, since the expectations, 
variances, etc. of its order statistics are simple and 
exact, and then to discuss the assumptions that are necessary 
before we generalize it for other distributions. 

a. Xen) for the Uniform Distribution 

het Xer-1)? Xery? X ere) be the order statistics 

from the uniform (0,1) distribution. In addition to the pth 
order statistic, a pooled estimator in the form of a linear 


combination of the (r-1), (r), (rt+l) order statistics will 


be considered. Here cur third candidate estimator is: 


Ze 





(x‘* puewsqd) 


; (4), ‘JoJeWTISY ezZeEptpuepg puoovesg “2g quYynold 


(T+4), (1) 


(Tt), 


X 


(* aTtzuend) 
OS9409 _ 
= —— Os) Zs 
fa 
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comer pe fe (e-t) Ce) A Gee) 


sil 


where Qt iS ey = es 


As mentioned earlier, if this estimator can be an unbiased 
estimator of xy? then a comparison of the variances of the 
ron order statistic (First candidate estimator), and this 

pooled estimator (Third candidate) will indicate which one 


has better performance. 


To maintain unbiasedness, we must have: 


E(X;,)] = EX yy] > oF 


es) 
= 
><il 
cn 
ry 
ae 
heel 
\ 


= Oo EX (p-1)! + B E[X(,] aay EX (y4y)]- 


If demand has a uniform distribution, then 


EX jl=¢ Gp+s GE)+y Gp 





This quantity is equal to E[X (7)! only when ao = y. 


= lS 
n+1 

One way to decide on the values of a, 8 and y is to choose 

values so that the variance is minimized. So, we would like 


to minimize 
Var (X (yy) = Var(a Xer-1) aro Xr) ne ey. Xtre1)? > 


Bubject to, 


i) 
j}—~ 


een? vy. and oa | 


By the constraints on a, 8 and y, these coefficients can be 


written in terms of a only, yielding 


Var (X/,y) = Var{a(Xp.4) * Xcpaqy) + (1-20) X(,y} 
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BY, 


Var (X(,y) = 0 “(Var (X )? + Var(X 


(r-1 (r+1)?? 


| oe Cov (KX (54? ¢r41)? 


+ (20-40%) (Cov (Xx y? + Cov(X 


(r-1)?*¢r (r)?* (r+1))? 


PRiGies er 4a“) Var (X (2) . 


Taking the derivative of Var(X/,y) with respect to a and setting 


it equal to zero gives, 


.. 2Var(X 7) - Cov yy Xr)? - Cov(X py oX rary) 


War (X ey_4)) * Var (hay) * ACOV ER 7 A (rer)? 7 
AL COW {Ky Xepy) * OVX ry Xpr4y) ~ VATA (yy )) 
(14) 
Having available the values of variance and covariances by 


equation (ll), equation (14) becomes 


Z2r(n-rtl) - (r-1)(n-rtl) - r(n-r) 


(r-1) (n-r+2) + 4r(n-rt+l1) + (r+1)(n-r) - 4(r-1)(n-rerl) - 
4r{n-r) + 2(r-1)(n-r) 


the n,r terms in the numerator and in the denominator cancel, 
resulting ina = 41/2. So, to minimize the variance of Xepy > 
the coefficients have to be a= y = 1/2, and B = 0. 

This says that when a population is uniform, a third 


reasonable estimator of xy 1S: 
Xr) = 1/2 (X 7) + Xre1y) : (15) 


a5 





yielding 


2r(n-r+l) - (ntl) 


Var (X/4) = ei ry 


If the ratio of the two candidate estimators is taken, 


<ileaal 9 : ér(n-r+1) ect te 
Var (X (yy) 2r(n-r+l) - (n+l) 


always. This implies that for a uniform distribution, 


Var (X ) > Var (Xpy)> and Xen) is the better estimator of 


(r) 
xX, In Figure 3, for n - 20, these variances are sketched. 
When r is close to bounds (either close to 1 or n), the 


difference between those two variances become larger. 


(Variance) 


0.015 


0.002 





Cm 


FIGURE 3. Variances of X and X o for Uniform(0,1) 
ASA nis {Gant 


Distribution,w y 
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These results apply, of course, only to uniform 
Penand distribution. In the following paragraphs the gener- 
alizations to any demand distribution will be discussed. 

b. Generalization of Xen) to any Demand Distribution 

For distributions other than uniform, it is 
difficult to claim that Xr)? as given by equation (15), is 
exactly unbiased. Also, the variance and covariances, as 
equation (13), will be functions of the density £,(X). Even 
though these two problems exist, some reasonable approxima- 
tions can be discussed. If we have both a large n and fairly 
flat densities, since eal (r), (rtl) are the closest 
neighbors among the other order statistic values, the densi- 
ties corresponding to these points would tend to equal each 
Other and during the calculation of qa, as in equation (14), 
would essentially cancel. 

If we have flat densities (where the three orlace 
are nearly linear), or if we have a large number of observa- 
tions (when the points are close together accordingly), 

Xen) would be approximately unbiased, too. Thus, it may be 
reasonable in some cases to use Xen) for other distributions 
besides the uniform. 

Incidentally, we would not like to try another 


pooled estimator of x, such as 


><Ill 


=O 


ce 2 cae gree = “Cry te ety 2 (reo) 


For this estimator, the assumptions which were made to permit 


Aa 





Xr) to-apply the distributions other than the uniform may 


no-longer be appropriate. For example, in equation (13), 

£ (Fy) (39)) and £)(F, > (Seq)) would probably have quite 
different values, because these are the densities of the two 
farthest order statistics, Xer-2) and Xr 42) in Ree 
Otherwise, if these densities were known, we would not use 

a nonparametric approach. 

So far, solutions of the Newsboy problem under 
"Uncertainty" and under "Risk" have been reviewed, and 
possible procedures for a transition phase between them 
discussed. Also in this chapter, three candidate estimators 
of the optimal inventory level So (or - were developed for 
use as alternative transition phase procedures. In the next 
chapter, the performance of these alternative transition 
procedures will be evaluated by a computer simulation of 
Newsboy problem decision making. Finally, from the simula- 
tion we will obtain information about what time to leave the 
minimax (Uncertainty) rule, for each transition phase pro- 


cedure. 
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IV. EVALUATION OF THE CANDIDATE ESTIMATORS 
OF OPTIMUM INVENTORY LEVEL, FOR THE 
TRANSITION PHASE OF THE NEWSBOY PROBLEM 
In Chapter III, the rationale for using a nonparametric 
approach for the transition phase of the Newsboy problem was 
discussed. It was shown that under Uncertainty the expected 
cost for a period was approximately linear in the mean square 
error of the optimum inventory level So (or x corresponding 


f° the pth 


quantile, and three candidate estimators of - or 
So were introduced. | 

In this chapter, the perfurmance of these estimators 
will be tested for several demand distributions by simulation. 
After explaining our reasons for using simulation, we will 
explain the structure of the simulation program used. In 
part C, the simulation parameters used in the program are 
described, together with the reasons for choosing these 
parameters. Part D illustrates the results of simulation. 
These results will be used later to compare the performance 
of the estimators of x and the transition phase procedures, 


and to see whether there were any apparent adverse effects 


of approximations or assumptions. 


A. RATIONALE FOR SIMULATION 

The purpose of the study was to find a transition phase 
procedure (i.e., which estimator of x, or Ds to use and when) 
during the transition from Uncertainty to conditions under 


Risk for the Newsboy problem, As explained before, the 
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performance of an alternative transition phase procedure 
depends directly on the performance of the candidate esti- 
mator of xy used, in that there is a one-to-one correspondence. 

A primary reason for simulation, of course, was to study 
the question of how long one should continue using the mini- 
max procedure (for Uncertainty) before changing to a transi- 
tion phase using a quantile estimate. Simulation over several 
demand distributions should provide insight regarding when 
(with respect to cost) this change in the Newsboy problem 
Should occur. Additionally, we wish to further examine the 
question of which of the candidate estimators to use. 

However, it was hoped that the estimators of x, or 35 
would come close to having desired order statistics proper- 
ties, being exact for the case of the uniform distribution 
(equations (9), (10), (11) and (12)). For example, for the 
th 


r order statistic Xr) it has been shown that 


Fy (EX 51) = 7 for a uniform distribution, but 


Fy (EX) 1) < 7 for convex, and Fy ELK (py 1) > 


A 


for 





bE 
n+] 
concave cumulative distribution functions of demand. [4]. 
Unfortunately, there is not an analytical measure to show 
how much the slack in these inequalities would change 
depending upon the amount of convexity or concavity of the 
demand distributions. 

Also, for the third candidate estimator of Ky? some 
additional assumptions were made that might not hold for 


every distribution. 
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In summary, our reasons for simulation are: 

1. To find the proper time (period) to begin the 
transition phase (to use one of the candidate estimators), 
leaving the "Uncertainty" case and the minimax decision rule, 

2. To provide information to supplement incomplete 
analytical knowledge of order statistics regarding popula- 
tion quantiles, 

3. To evaluate the impact of the approximations made 


in developing the third candidate estimator, Xery? over 


various demand distributions. 


B. DESIGN FOR SIMULATION 

The simulation program was written in the FORTRAN IV 
language. In every simulation run, inventory stocking 
decisions are made successively using one particular demand 


distribution and one particular value of the quantile 
C 


(P = oe ), say, demands generated from a normal parent 
Cc SO 
for quantile 0.3 ... etc. The demand distributions and 


various quantile values used in the simulation are given in 
the next section. Each run consisted of 50 successive 
decision periods starting with no demand information other 
than the range estimate, and replicated 60 times. 

As output values of the program, costs for each period 
were averaged over the number of replications. These average 
costs ner period were calculated for five different cases: 


1. The Uncertainty case, with the minimax decision rule 
being used, denoted by M. 


2. The ideal case (Risk), with the quantile value known, 
denoted by ID. 
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3. The transition phase case using estimator Kory? denoted 
DyouR 


4. The transition phase case using estimator Xory? denoted 
by R. 


The minimax (for Uncertainty) and ideal (for Risk) cases 
would give constant inventory levels, i.e., the same amounts 
all the way through the simulation run. For Uncertainty, the 
range of demand was approximately known and would not change 
unless it was revised, while for Risk eondions it was 
assumed that the distribution of demand, and thus the quan- 
tile, was known. 

Transition phase optimum inventory level values generated 
by aye Xen) and Xr) were subject to change from period to 
period until finally they converge to the value associated 
with the Ideal case. For the first two periods, Xery? Xery? 
Xen) would not lead to an inventory level, since they need 
at least two data points. 

The program was run for 50 periods, and replicated 60 
times. So, period costs for each transition phase procedure 
were averaged over 60 replications. Program steps can be 
ordered in the following way: 

1. Define the demand distribution yielding the quantile 
p, give the stockout cost Coo? the carrying cost Co» and the 
range of demand. 

2. Compute and set minimax and ideal case optimum 
inventory values. 

3. Generate a random demand value from the specified 


demand distribution. 
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4. Find the single-period cost for each case (minimax, 
Riny? Aen? Xr) and Ideal) by comparing their inventory 
_levels and the demand generated for this period. Record 
these costs. Find the next period inventory levels revising 
Xr)? Xr)? a) and go to the next period. 

5. Repeat steps 3 and 4 for 50 periods. 

6. Again, beginning from the first period, repeat 
step 5 for 60 times (replications), and at the end of 60 
replications find the average costs, and the sample variances 


of costs over 60 replications, for each period. 


7. Tabulate the results. 


The results obtained from this simulation will be examined 


in part D. 


C. SIMULATION PARAMETERS 

In this simulation a variety of possible forms of distri- 
butions, including continuous, discrete, symmetric, skewed 
(long-tailed), etc., were chosen. Demand distributions and 


their parameters used in this simulation were the following: 


Distribution Parameters 

Uniform the range over [0,30] 

Poisson with mean 10 units 

Chi-square with 10 degrees of freedom 

Lognormal generated by a normal parent 
having a mean of 3 and a variance 

of 1 
Normal with mean 35 and variance 100 
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For Da in the Minimax procedure, wp +3t1is used for 


C 
tailed distributions. Quantiles p = = » tested were 
Cc SO 


0.10, 0.30, 0.50, 0.70, 0.90, 0.95 for each distribution. 


ax 


These were determined so that Cus + C. was constant (in this 
simulation Coy + C. = 4.0 was used). So, for each {quantile, 
distribution} combination the differences of the expected 
period costs would come from the performance of the estimators 
of xy? as was shown earlier by equation (9). 


In the next section we shall present and discuss some 


of the results from the simulation. 


D. SIMULATION RESULTS 

In this part, simulation results will be evaluated and, 
for the Newsboy problem, a transition phase procedure will 
be proposed. In presenting these results we will identify 
the quantile value as p. Tables I and II give overall mean 
costs and standard deviations of cost. Then, Table III will 
show at what period the transition procedure should begin 
for a particular distribution, and quantile. Table IV will 
give relative efficiency of the best transition procedure 
to the Ideal case (decision under Risk). 

On the next page, Table I gives the overall mean values 
of the cost per period (over 50 periods and 60 replications) 


for each transition case, using C. + C. = 4,00. As can be 


O 


seen from Table I, there does not seem a significant differ- 


ence using R, R, or R, even though R looks quite consistent 


at the tail ( p= 0.95). One reason may be the fact that the 
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TABLE I 


Average Costs Over 50 Periods and 60 Replications 


having C,+C,,= 4.0. 


M=Minimax, R=X (yy > R=Xiyy> R=X (yy » ID=Ideal Cases 


Pp Uniform Poisson Chi-Square Lognormal Normal 


M 5.49 3.23 3.01 67.00 11.66 
R 6.35 Pes eel 14.50 8.09 
ile R 6.31 2.38 2.88 14.50 8.09 
R 6.30 2.38 OS 14.48 Se07 
ID 5.49 2.04 2.56 11.44 6.93 
M 12.68 4.60 5.92 255.96 20.04 
R 13.58 4.55 5.98 36.34 14.25 
0.3 R 13.55 4.52 5.96 36.23 14.18 
R 13.46 4.56 5.94 36.38 14.16 
ID 12.68 4.26 5.62 30.15 13.17 
M 15.08 Seon 11.08 339.86 se 7 
R 16.03 5.35 AG 51.22 15.80 
0.5 R 15.97 5.33 7.45 51.37 15.78 
R 15.91 5.34 7.45 51.35 15.73 
ID 15.08 4.98 6.98 43.49 14.80 
M 12.75 7.33 11.81 298.67 14.97 
R 13.75 4.74 etal 57.24 14.24 
wa7 OR 13.66 4.74 esa 57 ees 14.19 
R 13.60 A772 720 57.66 14.19 
ID 12.75 4.44 6.67 48.72 13.24 
M 5.39 4.37 6.26 132.42 9,33 
R 6.31 2.68 4.40 42.36 7.76 
0.9 R 6.30 2.67 4.38 42.38 7.76 
R 6.30 2.66 4.36 42.58 TTS 
ID 5.39 INS 3.71 35.91 6.70 
M 2.83 2.38 3.34 70.42 5.27 
R 3.77 els 3.00 31.80 5.04 
0.95 R 3.76 lens 2.99 31.74 5.02 
R 3.76 wes 3.00 32.44 5.04 
ID 2.83 1.41 I 25.47 3.91 
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4 
{ 
| 


difference between variance of Xr) and variance of Xepy> 

was not much (Figure 3). But if one wishes to see the dif- 
ference no matter how much it is, the following normalization 
by ideal costs may be considered: 


1. For a distribution from Table I, for every quantile 
R.-ID 
ik 
ID 
every R., where R:. is R, R or 
1 1 





p find , and sum this value over all quantiles p for 


‘all 


2. Repeat step 1 for every distribution. 
3. Pick the R;'s which have the smallest sum among 


Rm, &, R for every distribution as their best procedure. 


The results of this normalization were as follows: 


Uniform Poisson Chi-square Lognormal Normal 


Pee Ge ova 0.767 0.859 ess 0.839 
R 0.845 0.752 0.831 1.251 0.824 
R 0.828 1) HG 0.821 2G 0.816 


The best transition procedures were to use R for uniform, 
R for poisson, R for chi-square, R for lognormal, and R for 
normal distributions. To see the effects of the very right 
tail of a distribution on transition procedures, p = 0.99 
was used for normal and chi-square distributions. At this 
quantile the minimax costs were below R, R, and R costs. 
This says minimax can be used at the very tail of a distri- 
bution along the transition phase, providing very accurate 
prediction of the range of demand (say +5%). Otherwise, if 
we remember from Table I forp = 0.95, R, R and R were 


better than minimax approach. 


56 





From Table II (overall standard deviations over 50 periods 
and 60 replications), it is seen again that there is not much 
difference of standard deviations, even though R has a better 
trend (smaller sample variance of cost) among .R, R and R. 

From these two tables it can be said that there is not 
much difference of the costs yielded using R, R and R, but 
if an order is necessary it should be the one using R for p 
less than 0.95, and R at the tail (i.e., when p is greater 
than 0.95). 

During the simulation, as the number of periods increase 
the estimators converge to x, or the true So? while the 
minimax and ideal case inventory levels were constant. 

Since the values of X ery? Xr: Xr) usually change from 
period to period because they are updated by the most recent 
demand value, aha will generally be associated with a larger 
variance of cost than the minimix and ideal cases. This is 
reflected in Table II. Sometimes the minimax procedure had 

a smaller sample variance of cost (while possibly yielding 
the largest average cost) than the ideal case. This may be 
explained as follows: 

The ideal solution we are using relates to minimizing 
expected cost, not the variance of cost. These two minimi- 
zation points are generally different, depending on the 
distribution function. If the minimax solution is different 
from the ideal (expected cost) inventory level and closer to 
the minimum variance solution, it would yield a eoukier 


Variance of cost. Further derivations about the variance 
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TABLE II 


Sample Standard Deviations over 50 Periods 
and 60 Replications 


M=Minimax, R#X ry? R=X (1) > REX py» ID=Ideal Cases 


Pp Uniform Poisson Chi-Square Lognormal Normal 


M 3.15 1.26 ee nen 43.46 3.77 
R 6.36 2.69 Cae at) 27.46 9.89 
O.1 R 6.26 2.65 3.21 2a 9.70 
R 6.26 2.59 3.18 ea 9.53 
ID 3.15 1.62 2.26 15.90 6.62 
M 7.46 3.18 4.61 76.10 10.33 
R 9.00 3.60 4.93 57.64 2 alee 
0.3 R 8.84 3.53 4.89 57.69 Is 0) 4/ 
R 8.85 3.51 4.86 elk 11.99 
ID 7.46 3.16 4.42 45.93 10.95 
M 8.82 4.30 6.35 66.95 24 
R 10.31 4.14 6.37 84.29 12.73 
0.5. RK 10.23 4.12 6.33 83.91 12.68 
R 10.23 4.14 6.44 84.23 12275 
ID 8.82 3.84 5.80 PZ 12209 
M 7.30 3.47 4.87 46.79 9.62 
R 9.08 3.95 6.66 99.68 11.65 
Opa R 9.01 4.04 6.83 99.68 11.65 
R 8.87 4.09 6.76 99.23 oy 
ID (7.30 3.55 5.66 92.87 Od 
M 3.16 1.25 1.87 32.45 3.90 
R 6.59 3.28 5.72 99.99 9.10 
0.9 R 6.55 3.33 5.76 99.23 9.09 
R 6.55 3.29 5.69 98.05 9.01 
ID 3.19 2.57 4.10 94.42 6.34 
M 1.66 0.69 eee 31.66 27 
R 6.50 2.91 5.36 94.64 8.50 
0.95 R 6.50 2.95 5.34 92.85 8.47 
R 6.50 2.94 5.23 G27 4 8.46 
ID 1.66 1.49 325 84.97 Adal 
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minimization solutions of the Newsboy problem are given in 
Appendix B. 

The simulation results about the period to begin a 
transition phase are given in Table III. Shown are the 
first periods in which the average period cost (over 60 
replications) using the estimator was below that obtained 
using the minimax rule. In the Simulation it was observed 
that once costs using an estimator fell below the minimax 
costs, they stayed below for the remaining time periods. 
Simulation results showed that, on the average, at the fifth 
period we were supposed to begin a transition procedure, 
ceasing to use the minimax procedure. In other words, for 
the first five periods the large variability of candidate 
estimators of So (or x? caused a larger average cost than 
the minimax procedure yielded. As avery special case, if 
the distribution of demand is uniform, minimax and [Ideal 
case solutions would match each other and thus the uniform 


distribution is omitted from Table III. Also by coincidence, 


Cs 


i p = == is at or close to an intersection of the 
o=e SO 

implied uniform CDF of minimax procedure and the actual 
distribution function of demand, then the minimax procedure 
may yield the Ideal case solution. Otherwise, minimax in- 
ventory level has always a higher cost than the Ideal case. 
In Table III the large numbers are those specific values 
which relate to a quantile P that is close to minimax-- 


actual CDF intersection point. So it takes time for R, R 


and R to yield costs less than the minimax. No matter how 


Se 





TABLE III 


Periods to Begin Transition Phase 
Leaving the Minimax Procedure 


Quantile (p) Poisson Chi-Square Lognormal Normal 


Oa 6 6 5 6 
O55 Loe BO 4 4 
US 6 4 4 10* 
Dew 4 4 4 6 
oo 5 5 4 5 
0.95 6 5 4 6 


* Ideal-Minimax intercept points (or, close to these points). 
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long it will take to fall below the minimax cost (here, 
nearly at period ten) from simulation outputs it was observed 
that after the fifth or sixth period the cost difference 
between the minimax and the transition procedures was not 
large. This also can be predicted from the pattern values 

of Table III. The numbers in Table III are relevant for all 
candidate. estimators, i.e., R, R and R all fell below the 
minimax cost simultaneously at the same period. 

As a measure of efficiency for a transition procedure, 
the ratio of average cost induced by this procedure to Ideal 
average cost is used. Table IV gives these ratios for the 
average costs over second and fifth 10 periods for R (since 
R seems the promising transition procedure among the other 
ones). 

We will examine the data of Table IV for the second 10 
period, assuming that the minimax has been left and the 
transition phase has begun. We also believe that convergence 
toward the true value of Ideal inventory level So for the 
third, fourth, ... 10 periods will not be worse than the 
convergence accomplished during the second 10 period. 

If the transition-procedure (here, R) was operated at 
Micwtatile(@i.ec.5. 5 7 0-9), then the efficiency ratio tor the 
second 10 periods was on the average 1.16 and more (worse). 
As the quantile gets close to 1.0, the minimax procedure was 
yielding much better results than a transition procedure 


(assuming the range of demand is predicted quite well). 
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TABLE IV 


The Ratio of the Second and Fifth Ten-Period 

(ll to 20, 41 to 40) Average Costs to the 

Ideal Case (Under Risk) Costs, Using Xr) 
Transition Procedure 


Quantile (p) Uniform Poisson Chi-Square Lognormal Normal 


S 1.063 1.065 1.058 OSs lems 
re F 1.000 1.023 1.046 1.002 1.008 
S 1.065 1.059 1.050 1.040 1.066 
- 2 Og! Mesoca 1.003 1.004 1.008 
S 1.056 1.069 1.039 Loa = se onT 
Par EIS se 02 vols 1.004 1.007 
Ses, «1069 1.061 1.057 1.054 
si F 1.000 #£1.001 1.006 1.006 1.004 
3 Siig Sn ey 1.190 1.096 
= F 1.008 1.016 1.052 1.026 1.027 
eemieiom tos 1.211 1.240 1.183 
: 1.002 1.040 1.054 1.062 iOS 
Notes: . 


S stands for "SECOND," F stands for "FIFTH." 


When p is less than 0.9 overall ratio is 1.060 for second 
and 1.010 for fifth 10 periods. 
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For the rest of the quantiles where p is less than 0.9, 
this ratio for the transition procedure was almost constant 
1.060. To test the effects of quantiles and the effects of 
the form of demand distributions on the efficiency ratio 
(1.060), a two-way Analysis of Variance could be applied to 
the data of Table IV for quantiles less than 0.9. 

Different distributions were significant at a level 0.03, 
and the different quantiles were significant at a level of 


0.01. This result was explained that for C. + C. = 4.0 


O 
cost units, whatever the distribution of demand was and for 
any quantile less than 0.9, the ratio of the cost induced by 
transition phase procedure (Xp » here) to the Ideal cost 
(expected cost under Risk) would be practically constant at 
1.060 for second 10 period, and would not exceed this amount 
for the further periods (for example, at the fifth 10 period 
this ratio was 1.01). 

On the next part, as the summary of this chapter, and 
using the results obtained in it, a‘Single-Period inventory 
(Newsboy) problem transition phase procedure-is suggested. 
E. PROPOSED NEWSBOY TRANSITION PROCEDURE ALGORITHM 

(1) Approximate the bounds of Demand (Range). : 
(2) Set p= ee 
Ce tap 20. 90,760 80 Step so, otherwise go to step 4. 


(4) Use minimax for the first 5 or 6 periods. 


(5) For 6 < n < 50, use R (unless a reason for R, R), 
and go to step 10. 


Holt pes 0-95, go to step 8, otherwise go to step 7. 
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(7) 


(8) 


(9) 
(10) 


If step 1 is accurage, may use minimax for all periods 
and go to step 10, but if it is not accurate, go to 
step 8. 


Use minimax for the first 5 or 6 periods, then go to 
step 9. 


For 6 < n < 50, use R and go to 10. 


Analyze data. 


*May fit a probability distribution to the data and begin 
Decision Under Risk case. 


*May use R, R or R for n > 50 too. Since they converge 
quite fast, the risk of using them after 50 periods hope- 


X 
fHlly will be nearly equal to the risk of Ideal case. 
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V. CONCLUSIONS AND RECOMMENDATIONS 
FOR FURTHER STUDY 
For the Newsboy problem, a nonparametric decision pro- 
cedure that allows transition from Uncertainty to decision 


under Risk has been developed. Three estimators of the 


demand corresponding to the noe quantile were obtained, 
@ 

where p = C 22 Testing these estimators over several 
C SO 


distributions of demand with Newsboy cost criteria, it was 
found that excluding cases such as drastic skewness, pooling 
the (r-l)st and (rt+l)st order statistics for R provided 
essentially the best estimator, where R and R also worked 
well. 

It was also found that for the distributions tested and 
excluding Se values very close to 0.0 or 1.0, that 
the switching Newsboy decision making from minimax solutions 
to quantile estimator could be made, possibly at the sixth 
period. 

For further work, some parametric studies such as log- 
normal distribution, or beta are gamma families, depending 
upon whether the population of demand is finite or not, may 
be tried. Then, given ear the decision maker does not know 
what the distribution of demand is, which procedure, nonpara- 
metric or parametric, gives better response to Uncertainty 
may be compared. As another recommendation, some work can 


be afforded to improve nonparametric approaches further. 
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APPENDIX A 
ESTIMATION FOR POPULATION QUANTILES 
BY ,ORDER STATISTICS 
In this Appendix order statistics, especially those 
concerning population quantiles, will be briefly reviewed. 
For a further study of order statistics, references [1], 


fzieand [3] are cited. 
The Joint Distribution of n Order Statistics 


Order statistics theory is mainly based on continuous 
distribution functions. In reference [2] for discrete 
parents the joint distribution of order statistics is derived 
as well. 


Bet Xj 2Xg0--++9X, be a random sample from a probability 


distribution function FY. If Xejj2ials2,.-+-on is the , th 


largest value in this sample, then X < 


X 
(1) (2) 
are the order statistics of a random sample X 
th 


er eee Se € 


ee nO 


order statistic. Since Fy 1s 


a distribution function, by the probability integral trans- 


Here, X is Called the r 
(r) 


formation U= Fy (x) will have uniform distribution over [0,1], 
1.e., 
We EO a ec 


Bay (u) = 
0, otherwise. 


Letting Fy (x) = u with Fy (X) Seana Ome te Xs x, 

the following equations hold: (16) 
-1 -] 

Fy (u) = P(U<u) = P{F, (X) <u} = P[X<F, (u)] = FY [Fy (u)] =u 
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Since Fy (X) 1S a monotone increasing function and by 


equation (16) there is a one-to-one correspondence between 
Fy &qy) < Fy X(2y) < vec @ In) and Xe1) < X62) Ee eg 


Xin): Let fy be the density of By so that 


fy ee eo ee is the joint density of 
Ihe ave n 
@nessample X.5X,,..-,4.9. If each X-, 1=1,2Z,...,5n 
oar n 1 | 


1s independent and identically distributed, then 


i | 
}—! 


id (Cees Gh. it = Pee) 
X1oXor-004 9X i me n i 5 ete | 


Since there exists n! permutation of those sample values, 


Ih 
f CFs ye ee PCy (17) 
Xr XeayrreoXeny Y2 Yn i=17°X‘7i 


is the joint density of n order statistics. 


From the previous references cited: 


n! bg n-1 
fy & 290 ie @amienerm xo)! [1-Fy(y,)] fy Oy,)> 


(vr) 
(18) 


When Fy 1s a uniform distribution over (0,1), this equation 


turns to be, 


n! r-l n-r 
Rage = fealjnassyr bY four , (19) 


and the following sequence of equations are given, 


n! il 


BX ay] * Gerrmeryy fe Gu? au 


and 
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fer Grek )is.4(Et ln 


i) = 
E[X (+) es (n+k) (nt+k-1).... (+2) (nt1) 20) 
where k is a positive integer. For k=1: 
Re 
E(X yy] * get - (21) 


By equations (20) and (21) 


: Da ae Ze ar (n-r el) 
and, 


(r) ~ X¢sy! (n+1) (n+2) 





E[X 5 lo ae (23) 


From equations (6) and (8), 


Cov[X (1) »X(gy] = EX 1X5)! - EX, ]EIX(,)]> and 
= r(n-stl) 
Cov [X (1) »X 5)! Ok (24) 


Now, we wish a transition from the uniform distribution 


to a general distribution F Using the relations in 


x" 
equation (16), 


U(r) = F(X) or Xa) = Py Wrzy) 


th 


where Ver) is the r order statistic from Uniform (0,1). 


i 


7 Wen] % Fy* EW py) = Fy” Gey) 


Then, E[X x 


(r)! = E[F 


(25) 
Furthermore, the’ references given at the beginning of this 


appendix show that, for continuous distributions having 
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ie = inp,tl], and Osp, <P, -++ SP < 1 , 


the joint distribution of Seams) 2 (Gi) 222 Oe asymtotically 


tends to an %&-dimensional normal distribution with means 


Py» Poo--+> Py » and covariances, 
pi U-Py) oe 
coe Gu.) ro Ae x F(X , provided i < k, 
7 Pi PK. 
(26) 
So the variance of eae) will be, 
1 
p. (1-p.) 
Var [Xi J = = = 5 
i nt CCF 
Pi 
Having n, =T, it can be proven that, 
- Yn -X 
Xn.) XD. n (X(y) p,? 
U = = £85 (XX, ) N01) (27) 
ORE a)? yp. (1-p,) 1 


Equation (27) is quite important for this study. It 
says that, whatever the distribution function is, the pth 
order statistic, Xe? 1s approximately normally distributed 
around x (distribution free property). This also means 
that asymtotically Xr) 1s an unbiased estimator of x 


Let Gy (X) be the empirical distribution function of 


order statistics es 5) Sy 22 ony having, 


? ree cae, Gee X41) 
= Z ; 
Gy (X) =e if X05) < X < X (441) 
vay rs Xx 


(n) 
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In Reference [3] it is said, by the Glivenko-Centelli 
theorem that, Gy) converges uniformly to Fy (X) and nG p(X) 


is binomial (n, FY (X)), pee = ae 


pinGy(X)=j] = plGy(x) = 2) = (F)(FRQO]7 (1-Fy QI, 
(28) 


Point and Interval Estimation of Population Quantiles 


For point estimation of a population quantile, we assume 


that x is the ae quantile of the distribution function Fy, 
=1 

and P(X <X a x = eoG -f = X_. 

( p? See Je x ® : 

th 


Then, by equation (27) the r order statistic is a consistent 


and unbiased estimator of the pth quantile, where 
np , if np is integer 
_— Oso 
[npt1] , otherwise , and p = ~~ 
C2 eC 
so Cc 


For interval estimation of population quantiles, we are 


looking for r and s such that, 


<a. a, 4 


PUXery < Xp S A065)? 


“iV woere yY 1S percent Confidence. 
Now, 


ay, ime < 


Gay 5 ees) 
by equation (16) 


P(X = P(X) < X,) - P(Xpgy< X)> oF 


P(X ey) < x < Xesy) = PCFY (X25) <p) - PCF, (X55) Tee 


Since p = Py &) > by equation (28) 


Pi (Py pt(a-ppt 


S 


ita 


(7 )pta-pft- 
pw & al. 


] 
ios 


Xr) SS Se 


grt -i 
GO eG Dire = (29) 
1 
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To: summarize these derivations: 


Eh 


1. For large n the r order statistic is a consistent, 


h 


unbiased point estimator of the pt quantile for any distri- 


bile. on-. 


h 


2. Confidence intervals for the pt quantile as given 


by equation (29), distributed by the Binomial distribution 
(n,p). 

Sa es 1 gets larger, X er) converges to x, and also the 
empirical distribution function Gp (X) converges to actual 


cumulative distribution function Fy (X). 
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APPENDIX B 
EXPECTED COST AND VARIANCE MINIMIZATION POINTS 
OF THE NEWSBOY PROBLEM 

In this Appendix we will attempt to explain the differ- 
ence between minimization of the expected cost, and mini- 
mization of the variance of the cost, for the Newsboy 
problem. 

Decision makers' utility functions are involved with 
both expected cost (or return), and the variance of this 
cost. If a decision maker is risk prone, the variance of 
the cost may not affect him, and he just tries to minimize 
expected cost. On the other hand, if he is risk adverse, 
then he may give up some of the cost minimization value 
for a smaller variance (a trade-off between expected cost 


and variance of cost). As an example [5], let 
U(x) = x = bx4, be a utility function. Then, 
E[U(x)] = E[x] - bE[x“] = E[x] - b(Var(x) + [B(x)]%). 


If we try to optimize this value, we not only deal with 
expected value of x, but also variance of x. 

On the following pages, the variance minimization points 
of some distributions related with the Newsboy problem are 


derived. 


52 





Let Var(C) 2 aes - EC] - , and 


S 
E[C] = per (S-x) £, (x) dx + woe eo pe 
E cé = ae F S- de Only) be ae c 2 6 ¢x-s) 2¢ (n) dx 
[ory = C. " (S-x) fp = SO". D 


To minimize Var(C), with respect to 8: 


Dew ie 0, 1s necessary. 


dS 
a0 , 
Eee) = aE [C"] = 74310) cELC] = 0, (30) 
chen, 


[a @) 


2 S 
dE[C“} _ 2 : : fa 
as 2c : (Sex) 2 (G0) = AORe Lox S) £, (x) dx 


dE[C] | ° ie 
- 2E[C] aa = - 2E[C] (Cos a LOO ICRS C Sof se (Gslcbst 


- 2E[C]{C,Fp(S) - C,,(1-F)(8))} 





ds : 


S=5 


2 
since, el | = 2E[C] dE[C} 
S=S 
@) O 


then it yields 
2 S Dre 
C.tSF(S) - - £,(xjdx} - Coot fx Gs) Gee = Sip) a 


oO 


S 
(CS ICS) a i: Sse Ces : S(1-F,(S))}) 


(31) 
(Bf OCS. 8 Ge = G)- 





When this is applied to a uniform distribution over 


= 1 < < = 
edie. fy Ox) ye Xk Dee and FACS) D : 








(0, D 5 


nis max max 


then, the left hand side (LHS) of equation gives, 


D 
a ae 1 2. 2 “max 
LHS = S$ (C, - ©.)) 2D Secon Uso * 





Similarly, right hand side (RHS) of it yields, 


ORG s@ os oP max 
RHS = ((p=—59)s - C_.) (gs—— (C, + C,,) - C,8 + 22 BH). 
max max 


After we equate LHS - RHS and opened them, it was found that, 


50 Dax 
i) Scatter cerecemencnt 9 (32) 


: Co ss Ceo 


is the point that minimizes the variance of the cost for 
the uniform distribution. This point is exactly the same 
point that minimizes expected cost. 

Let us apply the same procedure to a beli-shape, say, 


Beta (a = 2, 8B = 2). Then 


r -1 -I 
a aeken °F q-xP += 6x(1-x) for 0 Sal 


i 


f(x) 


: 2 3 
F(x) = f 6u(l-u)du = 3X" - 2x”, and 
0 


ds 

L-Fp(x) = 1-3x?+2x9, fu £)(ujau = 4(1-4x9+ 3x4), 
x 

If uf) (u)du = 5 (4x° - 3x" ) 


Putting those values in equation (2), gives a 5 degree 


polynomial: 
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Ca 2 a4 Lae 
(Ge Gere, 5 mc Cae) S cmoaterae th Ge) one 
es 2 Z 2 2 3 Ze, 
: muce WO ‘ oC. C548 - (Gee i So i FECA, i Co) y 
(33) 


It may be difficult to solve this equation, but it is 


not difficult to see that it doesn't have the same solution 


& 
with SO 7 2 5 Loy a. Z 
55, 25, 5 OF 55, 25, p = Q, 


Cote So 


which is the minimum expected cost solution. 


In the exponential {i} case, 


=) X : a 
fy (x) = re d : Fy Ox) = l-e Kes 1-F yx) Sees and 


x - 
fu £(u)du = =[1-e eS eslas ei lee 
0 


oO 


fo (Ua 
De 


=[e**(1+ax)]. 


Again, inserting these values into equation (2) gives, 


Fy (S) - CA 
eae G4) 
c “so 
where FACS) = laa ee, 
C 56 ewe 
If we compare (34) with Fy(S,) = Gate l-e O 
C c “so 
or So = - 51n co): we will see that they are different 
c “so 


points. If we examine equation (34) further, it can be 


shown that S 7 So? 1.@€., the variance minimization point §S 


is closer to the tail than the minimum expected cost point So 


a 
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